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ABSTRACT 

Modeling molecular absorption by Earth's atmosphere is important for a 
wide range of astronomical observations, including broadband NIR photome- 
try and high-resolution NIR spectroscopy. Using a line-by-line radiative transfer 
approach, we calculate theoretical transmission spectra in the deep red optical 
(700 < A < 1050 nm) for Apache Point Observatory. In this region the spectrum 
is dominated by H2O, which is known to be highly variable in concentration on 
short timescales. We fit our telluric models to high-resolution observations of 
A stars and estimate the relative optical depth of H2O absorption under a wide 
range of observing conditions. We compare these optical depth estimates to si- 
multaneous measurements of Precipitable Water Vapor (PWV) based on data 
from a Global Positioning System (GPS) receiver located at Apache Point. We 
find that measured PWV correlates strongly with the scaling of II2O absorption 
lines in our spectra, indicating that GPS-based PWV measurements combined 
with atmospheric models may be a powerful tool for the real-time estimation 
of total molecular absorption in broad NIR bands. Using photometric measure- 
ments from the Sloan Digital Sky Survey (SDSS) DR8 database we demonstrate 
that PWV biases the calibrated r-z colors and z-band fluxes of mid-M stars but 
not mid-G stars. While this effect is small compared to other sources of noise 
in the SDSS z-band observations, future surveys like the Large Synoptic Survey 
Telescope aim for higher precision and will need to take time- variable molecular 
transmission into account for the global calibration of NIR measurements of ob- 
jects having strong spectral features at these wavelengths. Empirical calibrations 
based on PWV may be immediately applicable to ongoing efforts to make mmag 
differential measurements of M stars to detect transiting exoplanets. 



Subject headings: Astronomical Techniques 
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1. Introduction 

For many astronomical measurements it is necessary to understand and precisely mea- 
sure the attenuation of li ght by Earth's atmosph ere. While attempts have been made to do 



this in an absolute sense (iHayes fc Latham! Il975l ). in general astronomical measurements at 



optical and NIR wavelengths are fundamentally differential, measuring flux relative to either 
a standard star or comparison stars within a given field of view. As shown in Figure (H 
across the optical and infrared Rayleigh scattering by molecules, Mie scattering by aerosols, 
and line absorption by molecules leads to atmospheric tr ansmission that is a rapidly varying 
function of wavelength and can also vary rapidly in time (jstubbs et al. 2007 ). Massive imag- 



ing surveys like the Large Synoptic Survey Telescop(E(LSST), Dark Energ y Surve}i|(DES), 
Pan-STARRs|f|, and Hyper Suprime-Cano seek to produce broadband photometric measure- 
ments across the UV and near infrared (NIR) that are stable and uniform to better than 
1% over periods of years. Even if the response of the telescope to monochromatic light o f 



a given intensity is known very precisely (e.g. IStubbs fc Tonryl |2006| . IStubbs et al.ll2010l ). 
global calibration of photometric measurements requires correcting flux for atmospheric vari- 
ations. Similarly, precise radial velocity measurements at NIR wavelengths require a detailed 
understanding of atmospheric transmission so that subtle Doppler shifts in stellar spectral 
features can be separated from changes in telluric absorption features. 

The flux in ADU measured by a single-element detector from an astronomical source, 
r, can be expressed as a function of the wavelength-dependent properties of the telescope, 
atmosphere, and source as 

poo 

F = AAT / FxT{\)S{\)dX (1) 



JO 

where A is the effective collecting area of the telescope, AT is the exposure time, T(A) is the 
transmission of the atmosphere as a function of wavelength, S{X) is the response function of 
the telescope in units of ADU/photon (including the effects of optics, filters, and detectors), 
and Fx is the flux density of the source in units of photons s^^cm^^cm^^. Historically, for 
broadband photometry the band-integrated value of the product S{\)T{\) is measured for 
a given night using a combination of observations of a uniformly illuminated screen (the 
"fiat field") and observations of photometric "standard" stars. Using this methodology it 
is possible, under excellent observing conditions, to measure the absolute brightnesses of 



^\protec-t|ittp://www. lsst.org/ 



^Htp: / /www. darkenergysurvey.org/ 
^ttp://pan-starrs. ifa.hawaii.edu/public 
ihttp://www. naoj.org/Projects/HSC/index. html 
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objects with a precision of 1%, though these measurements are only as accurate as the 
measurements of the standard stars themselves. 

Temporal changes in T(A) and >S'(A) become important when the target and standard 
are not observed simultaneously. For many applications, such as precise photometry of 
transiting extrasolar planet hosts, only the relative brightness of an object, as compared 
to other nearby stars measured simultaneously, is required. This di fferential photome tric 
approach results in such high photometric precision, often 0.1% (e.g. IChan et al.ll201ll ). in 
part because the temporal variations in T(A) and 5'(A) are intrinsically removed. Given 
multiple observations of a sufficiently large number of objects, it is possible to generalize this 
differential approach and apply it to data from large-area survey s to internally calibrate the 
brightnesses (and colors) of objects. iPadmanabhan et al.l (|2008[ ) applied a technique of this 
type called Ubercal to the Sloan Digital Sky Survey (SDSS) photometric measurements and 
demonstrated that 1% global phot ometric precision can be achieved for a survey of millions 
of objects spanning years (see also llvezic et al.ll2007l ). 



Differential calibration strategies may not account for the fact that Fx varies strongly 
with stellar spectral type (or quasar redshift), resulting in higher-order effects that can arise 
when measuring relative fluxes between objects having very different Spectral Energy Distri- 
butions (SEDs). This is particularly important at NIR wavelengths, where strong molecular 
bands due to H2 O, CH4, and CO2 dorainate T(X) and Fa niay have strong abso rption or 



emission features (IBailer- Jones &: Lamnj|2003l : iBlake et al 



2008 



Mann et al.ll2011f ). In these 



cases, temporal variations in T(A) are not effectively removed through differential measure- 
ment, and residuals at the level of 1 mmg are expected in differential photometry of cool 
stars. Ultimately, the highest precision will be achieved for a survey like LSST with the help 
of detailed models of T{\) that are contemp oraneous with broad band photometric observa- 
tions and localized to particular sight lines. iBurke et al.l ( I2OIOI ) demonstrate an ambitious 
strategy whereby simultaneous spectrophotometry of calibration stars in LSST flelds can be 
combined with radiative transfer calculations for the atmosphere above the observatory to 
measure all important components of T(A). 

From the observational point of view, the challenges of precise radial velocities at NIR 
wavelengths and the global calibration of NIR photometry in large surveys are closely related 
since both require contemporaneous modeling of transmission through Earth's atmosphere. 
Precise stellar radial velocity measurements at wavelengths longer than 600 nm require a 
detailed understan ding of the impact of absorption due to molecules in Earth's atmosphere 
( IWang fc Gell201ll ). These telluric absorption features themselves can also be used as a 
simultaneous wavel ength referen c e, a m easurement technique co nceptually similar to th e 
I2 cell described by iButler et al.l (Il996l ). though flrst proposed by iGriffin fc GriffinI ( 1l973l ). 
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The use of telluric lines as a wavelength reference has a long history in the literature (e.g. 
Smitlj|l982l . iDeming et al.lll987l iHatzes &: Cochranlll993l ). and recently several authors have 
den ionstrated that t e lluric li nes are a viable v el ocity reference for pla net se arches, includ- 



Seifahrt fc K^ui feoosh . Isiake et al.l toid ). iFigueira et al.l and iFigueira et al.l 



This technique requires an excellent r nodel of T(A), whi ch can be derived empirically 



from high-resolution observations of the Sun (IBlake et al.ll201Cll ) or calculated directly based 
on as sumptions about the composition and structure of Earth's atmosphere (ISeifahrt et al. 
2010h . 



Atmospheric constituents like CH4, CO2, and O2 are well-mixed from Earth's surface 
through the mesosphere and have seasonal variations in total column density. For wave- 
lengths 900 < A < 2500 nm, absorption by H2O becomes very important, having large 
optical depth at some wavelengths. Modeling these absorption features poses a significant 
challenge since the total water column can have large and rapid variations. We demonstrate 
that Precipitable Water Vapor (PWV) estimates derived from a Global Positioning Sys- 
tem (GPS) receiver can be used to calculate theoretical transmission spectra having no free 
parameters and resulting in excellent fits to astronomical spectra in the wavelength range 
900 < A < 1000 nm. In a region of non-saturated H2O lines we find fit residuals < 2%, likely 
dominated by residual fringing signals in our spectra. We find that GPS-based PWV esti- 
mates are also useful for estimating total transmission in the SDSS z band and demonstrate 
that changes in PWV bias photometric measurements of mid-M stars and determine an 
empirical correction factor for the z-band photometric measurements in the SDSS database. 



2. ARC Echelle Observations 

As part of an ongoing program to measure the radial velocities of M stars, we obser ved 



74 A star telluric standards with the ARCES echelle spectrograph (IWang et al.l l2003l ) on 
the 3.5 m telescope at Apache Point Observatory (APO). A single ARCES observation 
contains data in more than 100 spectral orders spanning 360 nm to 1010 nm at a resolution 
of R~30,000. These A star observations were collected during 2010 under a wide range of 
observing conditions and at airmass (AM) up to AM=2.0. Exposure times were tailored 
so that S/N per pixel in the extracted spectra would be ~ 500 at A ~ 800 nm. The 
ARCES spectral orders are closely spaced and marginally sampled in the spatial direction, 
requiring special care in the spectral extraction process so as not to induce aliasing. We 
collected a large number of bias and quartz lamp frames on each night and used the quartz 
lamp frames to define the positions of the centers of each of 75 spectral orders, which run 
approximately parallel to the detector rows. After subtracting an average nightly bias, we 



- 5 - 



estimated scattered light in both the stellar spectra and quartz lamp spectra by fitting a low- 
order polynomial surface to the 74 inter-order minima of each column. Next, we extracted the 
stellar and quartz lamp spectra using an aperture extraction algorithm with fractional pixel 
weighting at the edges of the aperture. Since the orders are so narrow, it was necessary to 
fiat field the ARCES data using the ID extracted quartz lamp spectra, so we made a master 
flat for each spectral order for each observing night by averaging the individually normalized 
ID quartz spectra. For the red orders of interest here there is large amplitude fringing in 
the ARCES spectra, which is partially removed from the stellar spectra by dividing by the 
normalized master flats. There are slight (~1 pixel) offsets between the positions of the 
spectral orders in the quartz and stellar frames, resulting in incomplete fringe removal. 



3. Atmospheric Models 

Modeling telluric H^Q absor ption in our A star spectra requires high-resolution telluric 
templates. In iBlake et al. an observed atmosph eric transmission spectrum, derived 

from solar observations by iLivingston fc Wallace! (Il99l[ ). was scaled to fit telluric absorption 
features due to CH4 around 2300 nm. We found this approach to be problematic for absorp- 
tion by H2O between 800 nm and 1000 nna, so instead we created theoretical transmission 
spectra. Closely following iGordley et al.l (119941 ) we developed a custom line- by-line radia- 
tive transfer code to calculate template spectra that contain only molecular absorption (no 
scattering, no line mixing, no H2O continuum absorption) gi ven an atmospheric m odel and 
molecular transition data from the HITRAN 2008 database (IRothman et al.l 120091 ). Follow- 
ing the Beer-Lambert law for extinction through a medium by a single transition of a single 
absorber, the transmission, T^, of the Earth's atmosphere at frequency u is given by 



T 



(2) 



where cr(z/) is the frequency- dependent cross section and u is the mass path 



u 



qPx 



(3) 



where x is the path length through the atmosphere, P is pressure, T is temperature, k is 
Boltzman's constant, and q is the Volume Mixing Ration (VMR), which is defined as the 
fraction of the particles in a volume that are the given constituent. Earth's atmosphere 
contains many gasses and the properties of the atmosphere, including P and T, change with 
altitude and time. The transmission through the atmosphere can be calculated by breaking 
the atmosphere up into a number of levels in the vertical direction, z, within which the bulk 
properties are assumed to be constant. The total transmission due to a number of different 
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V{iy) = — — / 5 ^dt (5) 

JJ 7/ 



absorbing gasses, (7, each with a number of absorption hnes, z, can then be calculated 

The shapes and strengths of the absorption lines depend on the cross section, o"(z/). Ab- 
sorption lines in Earth's atmosphere have intrinsic shapes both due to the bulk properties of 
the gas mixture and the quantum mechanical properties of the transition. The approximate 
shape of an absorption line is given by the Voigt profile, Vlu), which is the convolution of line 
broadening kernels due to Doppler and pressure (Lorentzian) broadening mechanisms. At 
low pressure the motions of the absorbing molecules dominate the broadening and at higher 
pressure collisions between molecules come to dominate. Both regimes occur at different 
levels in Earth's atmosphere. The Voigt profile is defined 

.0 

aD7r3/2 aD J [{u - z/J /an - t]' + [aL/a^] 
where an is the Doppler half width of the line, is the Lorentzian half width of the line, and 
Uc is the line center, each of which can depend on atmospheric temperature and/or pressure. 
The cross section cr(z/) is a function of the line broadening as well as the line intensity, 
S = S^Q, where Q is a scaling factor based on the partition function of the molecule and the 
temperature of the gas. The total transmission can now be expressed as Tj, = e"^" where t^, 
is 

= 5^ 5^ Ug,Z Vg,i,,{u)Sg^i^, (6) 

z g i 

We retrieved line parameters including half-widths, strengths, line centers, and pressure 
shifts of those line centers, for transitions of H2O in the wavelength range 600 to 11000 nm 
from HITRAN 2008 and calculated atmospheric transmission at zenith in wavelength steps 
of O.OlA across a plane parallel atmosphere with 30 layers logarithmically spaced between 
km and 80 km. 

This radiative transfer calculation requires a model of the structure and composition 
of Earth's atmosphere, including temperature, pressure, and chemical composition, as a 
function of height. We used the NASA- MSFC Earth Global Reference Atmospheric Model 



(Earth-GRAM 2010; I Justus et al.ll2004j and references therein) to approximate the average 
properties of the atmosphere above APO, located in southeastern New Mexico at an altitude 
of 2,780 m. Earth-GRAM combines a large number of historic atmospheric measurements 
with modern meteorological models to predict atmospheric properties, and the statistical 
distributions of those properties, as a function of altitude, location on the Earth, and month. 
While Earth-GRAM is often used for engineering applications, such as estimating the effect 
of Earth's atmosphere on trajectories, it is also a powerful resource for quantifying the impact 
of the atmosphere on astronomical observations. 



- 7- 



4. Spectral Fitting 



Following iBlake et al.l (|2010[ ) we adopt a forward modeling approach to fitting our 
ARCES A star observations. We focus on order 3, spanning the wavelength range 970 < A < 
986 nm, where there are numerous unsaturated H2O lines and the A stars are intrinsically 
devoid of spectral features. The model, B{i), begins with the atmospheric transmission spec- 
trum based on the average Earth-GRAM model, Tq, scaled for water optical depth, r, which 
is a function of zenith angle (or airmass) and water content of the atmosphere. This scaled 
spectrum Tq is convolved with an estimate of the spectrograph hne spread function (LSF), 
interpolated onto a lower resolution pixel grid that is defined through a polynomial relation 
between pixel, i, and wavelength, W{i), and finally normalized by a continuum function, 
N{i), that is a polynomial in i 



B(i)= [{TS{X)(^LSF)\wir)] x N{t) 



(7) 



where ® indicates convolution and I indicates "evaluated at" . This model has 10 free param- 



eters : four for N{i), three for W{i), two for the LSF, which is a ID Moffat function ( iMoffat 
19691 ). and r. We fit this model to order 3 of e ach A star spectrum using an implemen ta- 
tion of the AMOEBA downhill simplex method jNelder fc Meadelll965l . IPress et al.lll992[ ) to 



minimize x^- Examples of fits are shown in Figures [2] and [31 The agreement between our the- 
oretical telluric template and the observed spectra is excellent overall, with typical residuals 
of 1 — 2% and no evidence for degradation of the fits in the centers of prominent absorption 
lines for regions where the line density is modest and the average line optical depth is < 1. 
As shown in Figure [31 in regions having numerous lines with large optical depths the overall 
quality of the fits deteriorates and the models exhibit significant discrepancies in the cores 
of some absorption features. In all cases the fit residuals do exhibit slowly varying features 
that we attribute to imperfect correction for CCD fringing, which has a large amplitude at 
these red wavelengths. Uncertainty in the HITRAN H2O line parameters may also lead to 
discrepancies between our models a nd observations and we note that incorporating empirical 
line parameters, such as those from lAlekseeva et al.l (l2010bl ). may help to improve the overall 
quality of the fits. 



5. PWV Measurements 

We found that overall the best-fit model parameters had very small scatter, demonstrat- 
ing that ARCES is physically stable, with small shifts in wavelength solutions over the period 
of a year. However, we expect that r, which is presumably related to zenith a ngle and the 



amount of H2O in the atmosphere, should vary widely since H2O, unlike CH4 (iBlake et al. 
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20101 ). is in- homogeneously distributed in the lower atmosphere and its concentration can 
change by more than 10% in under an hour. As a result, our best fit values of r spanned 
the ranged 0.1 < r < 0.9. The amount of H2O in the atmosphere is related to Precipitable 
Water Vapor (PWV), a measurable quantity that should be proportional to the optical 
depth to absorption by this molecule. Historically, PWV, which is critically important for 
mm and sub-micron observations, has been measured with balloon-borne radiosondes or dif- 
ferential radiometers operating in the mid-IR. In the 1990s it was demonstrated that data 
from a multi-band geodetic quality Global Positioning System (GPS) r e ceiver and a high 



accuracy baromete r can be combined to derive PWV (iBevis et al.l Il992l . iBevis et al.l Il994 



Rocken et al.lll995l ). making these devices extremely important for studies of the time and 



spatial dependence of water vapor for atmospheric science applications. 

The index of refraction of the atmosphere, and therefore the speed of propagation for the 
microwaves used by GPS, is a strong function of pressure, temperature, and water content. 
Both the ionosphere and neutral atmosphere induce delays in signals from GPS satellites 
that are measured precisely by the receiver stations. The ionospheric delays are frequency 
dependent, so GPS signals at multiple frequencies (typically 1.2 and 1.6 GHz) can be used to 
correct for this effect. The total delay at zenith (TD) induced by the neutral atmosphere can 
be broken into two components: "zenith hydrostatic delay" (ZHD), which is due primarily 
to dry air, and "zenith wet delay" (ZWD), which is due to H2O vapor. The ZHD is a slowly 
varying function of time and can be calculated given a barometric pressure measurement to a 
precision better than 0.01%. A measurement of the TD and an estimate of the ZHD therefore 
provide an estimate of the ZWD (ZWD=TD-ZHD), which is proportional to PWV with a 
constant of proportionality that is a function of surface temperature. While it is possible to 
calculate this constant of proportionality to make absolute measurements of PWV, relative 
measurements of PWV can be made with quoted precision better than ~ 0.2 mm given 
contemporaneous surface temperature and barometric pressure measurements. 

A large number of GPS-based PWV monitors are maintained around the world. Data 
from several receiver network, irr North America are maintairred by UCAR through SncnuNeg 
including a station operated by UNAVCO located in close proximity to the APO 3.5 m tele- 
scope. Using the SuomiNet web interface we downloaded PWV estimates at 30-minute 
intervals from the APO station for the 2009 and 2010 seasons and interpolated the PWV es- 
timates to the times of our A star observations, excluding cases where there was not a PWV 
estimate within 1 hour of the time when an A star spectrum was obtained. Using multiple 
linear regression, we fit for our measured values of r, the scale factor specific to our telluric 



' ]\prot eclj^ttp: / / www.suominet.ucar.edu/ 
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template, as a function of PWV and (AM-1) for PWV< 8 mm. We found that measured 
PWV and AM are strongly correlated with r and determined a best-fit linear relation 

r = 0.036 + 0.098 x PWV{mm) + 0.36 x {AM - 1) ± 0.055 (8) 

The results of this fit are shown in F igure HI W e note that this resul t is qu alitatively similar 
to those reported by iQuerel et al.l (120081 ) and iThomas-Osip et al.l (120071 ) who found that 
PWV estimates derived from radiometer observations strongly correlated with the measured 
depths of H2O features in optical spectra. As shown in Figure |5l our A star observations 
spanned a relatively narrow range of PWV compared to the overall distribution of PWV 
at APO, but we are biased toward lower PWV since periods of large PWV are likely to be 
unusable for astronomical observations. In practice, the coefficients of a relation like the 
one given in Equation 8, which is based on telluric models calculated using an "average" 
atmospheric profile for APO, will need to be calibrated once for a given observing site. Once 
this is done with a set of A star observations, the technique presented here should be broadly 
applicable to NIR observations at sites other than APO. 



SDSS Observations 



Our fits to the A star spectra demonstrate that given an estimate of PWV we can gen- 
erate a theoretical model for atmospheric transmission in z band that has no free parameters 
and matches the observed telluric absorption around 980 nm to better than 2% on a per- 
point basis. Integrated across a wide observing band, the discrepancy in total transmission 
between our models and the A star observations is significantly smaller, less than 0.2%. This 
opens up the possibility that GPS-based PWV estimates can be used as input to theoretical 
radiative transfer calculations to enable the precise estimation of molecular absorption inte- 
grated across NIR bands in real time. When combined with information about the spectra 
of the objects being observed, this technique could potentially lead to significant improve- 
ments in global NIR photometric precision fo r surveys like LSST and augment an overa ll 
calibration strategy like the ones outlined by iBurke et al.l (120 lOf ) and IJones et al.l ( I2OIOI ). 
We test this hypothesis by considering z-band photometric measurements of M stars in the 
SDSS photometric database. 



The SPS S (|York et al.l 120001) made an imaging survey of the sky using a la rge-format 
CCD camera ( iGunn et al. I ll998h mounted on the Sloan Foundation 2.5 m teles cope (iGunn et al 
2006^ at APO to image the sky in five optical bands - u, g, r, i, and z (IFukugita et al.l 



19961) . The imaging data were reduced by a set of softwar e pipelines that pro d uced a cat 



alog of objects with calibrated magnitudes and positions (ILupton et al.l l2002l iHogg et al. 
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200ll . IPadmanabhan et alJboOSf ). SDSS Data Release 8 (DR8; lAihara et al.lboilh includes 
photometric measurements in five bands of over 260 million stars over an area of more than 
15,000 square degrees. The GPS-based PWV monitor at APO (where both the SDSS and 
3.5 m telescopes are located) was installed in 2007 and began regular operations in 2008. 
Unfortunately, nearly all the SDSS imaging was acquired prior to this. A similar GPS-based 
PWV monitor operated by NOAA§ is located at White Sands (WS), New Mexico, 50 km 
southwest of APO at an elevation of 1200 m (1580 m below APO). We interpolated the PWV 
estimates from WS during 2009 to the times of those from APO having PWV< 10 mm, only 
considering APO measurements within 30 minutes of a WS measurement. We found that 
the APO and WS PWV measurements are highly correlated (Pearson's r=0.85), though 
with significant deviations from linearity at PWV< 1 mm. We adopt an empirical relation 
between WS and APO by interpolating the average PWV data given in Table [H While the 
scatter in this relation is large, particularly at low PWV, it is sufficient for considering the 
relationship between the overall statistical properties of SDSS z-band photometry and PWV. 

The photometric measu rements in the SDSS pR8 d atabase have undergone the calibra- 
tion procedure described in IPadmanabhan et al.l (120081 ). which has been shown to achieve 
global photometric precision of ~ 2% in z band, and ~ 1% in g, r, and i bands. In z band, in 
addition to the components of atmospheric extinction that vary smoothly with wavelength, 
such as Rayleigh scattering by molecules, Mie scattering by aerosols, and "gray" absorption 
by thin clouds, line absorption by H2O becomes important. This means that the extinction 
for a given object depends on the SED of the object, which is not the case for the other 
bands. Ubercal is fundamentally a differential approach, where flux is calibrated relative 
to a large ensemble of objects having a wide range of SEDs that, when averaged, at NIR 
wavelengths have a smooth SED like a solar type star. Given this, we expect systematic 
residuals in the SDSS z-band measurements for objects with complex NIR SEDs, such as 
M stars, that will correlate with PWV. We selected two samples of M stars from the SDSS 
DR8 databas^ to investigate this effect. 



The first s a mple is based on the spectroscopically-classified catalog of M stars produced 
by I West et al.l (1201 ll ). We selected 7,300 objects having spectral types from M7 to M4 
and brighter than i=17.0 and retrieved the CLEAN photometric measurements of these 
objects from the STAR view in the DR8 context. We estimated PWV at the time of each 
observation based on the WS data and rejected epochs where no WS data was available. 
We also rejected outliers in (r-i) vs.(r-z) color space by excluding objects with (r-z) colors 



IXprot eel ^ttp://www. gpsmet.noaa.gov/ 



1\protec-t[ittp:// skyserver .sdssS . org /dr8/en/tools/ search / iqs . asp 
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falling more than 0.2 mag from a linear fit to the stellar locus and calculated A(r — z), the 
color of each object relative to the stellar locus for a given (r-i). We also calculated the 
mean, with outlier rejection, of A(r — z) in each of the six camcols and found significant 
offsets which we subtracted off. The relation between A(r — z) and PWV is shown in 
Figure [6l The best fit linear relation between change in (r-z) color and PWV is A (r — 2;) = 
(-0.0019 ± 0.0002) X {PWV - 3.3 mm) mag. 

The second sample contains objects in Stripe 82, a 300 square degree area along DEC=0 



that was repeatedly observed by SDSS (llvezic et al 



Ivezic et al. 


2OO7I: 


Bramich et al. 


2008 



20081 ). Based on the stellar locus defined in lCovey et al.l (120071 ) we selected a r andom 3,000 



G5-F5 stars and 2,500 M7-M4 stars from the Stripe 82 standard star catalog of llvezic et al. 



( 120071 ). restricting our samples to objects brighter than z=17.0. While parts of Stripe 82 
have been observed up to 80 times, only a portion of those visits are included in DR8 
( lAihara et al.l 1201 if ). While all the visits are included in DR7, these data were not globally 
calibrated using Ubercal. We selected photometric observations of our Stripe 82 targets 
from the DR8 PHOTOOBJ view requiring all the photometric fiags required for CLEAN 
status except for PRIMARY. We culled observations at epochs where no WS PWV data 
was available, as well as observations wit h colors offset more than 0.1 mag from the (g-i) 
stellar locus defined in ICovey et al.l ( 120071 ). In cases where observations of the same object 
spanned multiple camcols, we chose a single camcol for that object and rejected observations 
from other camcols. We also excluded objects with fewer than five remaining observations 
at PWV< 4 mm. The final samples contain 14,597 observations of 983 G5-F5 stars and 
25,005 observations of 1,553 M7-M4 stars. For each object in each sample we calculate (z)^, 



the average of the z-band magnitudes only for PWV< 4 mm, and then Az 



(2;) 4. The 



relation between Az and PWV for both samples is shown in Figure [71 As expected, the effect 
of differential extinction due to H2O is small for G and F stars, which have SEDs similar 
to the SED of the average point source used in the photometric calibration, but the effect 
is large for M stars, with Azh^q = (0.0012 ± 0.00008) x {PWV - 2.6 mm) mag. The slope 
of this relation is consistent with that of the A(r — 2;)-PWV relation, indicating that our 
results obtained using two different sets of observations are very similar. Considering that 
the distribution of PWV for the Stripe 82 observations is roughly uniform from 1 mm to 10 
mm, this is potentially an important component of the global calibration of SDSS M star 
observations. However, taking the standard deviation of the empirical corrections for all the 
SDSS measurements indicates that this is a small effect overall: {AZH20) = 0.003 mag. For 
the M star s we recalculated Az after applying this empirical PWV correction. Using the 
F-statistic (jPress et al.lll992l ) we compared the variances of Az before and after the PWV 
correction and found no statistically significant reduction in the variance of Az {P{(j\ = 
al) = 0.44). 



- 12 - 



Using the Stripe 82 sample we can also investigate changes in overall z-band trans- 
mission with PWV by considering uncalibrated flux measurements. Using the PSFFLUX 
and NMGYPERCOUNT entries in the PHTOTOOBJ view we can convert the calibrated 
fluxes retrieved from the DR8 database back into uncalibrated, raw fluxes. Using the same 
procedure as for the calibrated Stripe 82 data, for each object in the G5-F5 sample we 
calculate C^z^aw = Zraw/ {zraw)^ and Avraw = fraw/ (rraw) 4- As a rcsult of cloud covcr and 
variations in the optical depth of molecular scattering, there are variations at the level 
of ~ 10% in the uncalibrated r and z fluxes, but Avraw and Azraw are highly correlated. 
We flt a linear relation, F{Arraw), for Azraw as a function of Avraw and then calculate 
a' Zraw = Azraw — F{Azraw)- Siucc the SDSS r band is free from strong H2O line absorption 
and probes instead the transmission variations due to scattering, we use A! Zraw as a proxy 
for the molecular component of absorption in z band. The change in raw z flux as a function 
of PWV is shown in Figure [3 The best flt linear relation between molecular absorption and 
PWV is /^Zraw{%) = (-0.0096 ± 0.0001) x {PWV - 4.3 mm). 

We compared these empirical relations between z raw flux, z magnitude, and the (r- 
z) colors of M stars to theoretical estimates of the impact of changing PWV on SDSS 
z-band transmission . We numerically integrated the synthetic stellar spectra described by 



Brott fc Hauschildtl ( 120051 ) across the SDSS passbandqfl including our H2O absorption tem- 



plate for APO scaled following our derived t{PWV) relation given in Equation 8 and assum- 
ing AM=1.2. We focused on stellar templates representative of a mid-M dwarf (Tcfr=3000 
K, log(^)=5.0) and a mid-G dwarf (Teff=6000 K, log(^)=4.5). For values of PWV in the 
range 0.5 to 12 we calculated the expected z-band and r-band fluxes for each template as 
a function of PWV, normalizing to the expected flux at PWV=3 mm. While the expected 
z-band fluxes change by more than 8% when PWV increases from 1 mm to 10 mm, to flrst 
approximation these variations are removed by the SDSS photometric calibration. If the 
appropriately weighted average of the SEDs of all SDSS stars is well-approximated by a 
G-star spectrum, then the differential fluxes we calculate in r and z bands should match 
the variations we measure in SDSS z-band photometry as a function of PWV. These the- 
oretical expectations are shown as dashed lines in Figures |6] and [7] and are in very good 
agreement with our measurements. We experimented with including different Aerosol and 
Rayleigh scattering models in these calculations, but as functions that vary smoothly with 
wavelength they had a negligible effect on the predicted differential M dwarf colors, fluxes, 
and differential magnitudes. 



^Xprotecl] 



^ttp: / / www.sdss3.org/instruments / camera. php 
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Discussion 



Empirical corrections for extinction due to H2O absorption may be important for a 
wide range of spectroscopic and photometric astronomical observations at deep red and in- 
frared wavelengths. In many cases, contemporaneous PWV measurements derived from GPS 
may provide an important external parameter that can be used to correct for these effects. 
Precise differential m easurements of ste l lar brightness can be used to find rocky planets tran- 



siting low-mass stars ( iBlake et al.ll2008t iNutzman fc Charbonneaull2008l : ICharbonneau et al. 



20091) as well as measure the properties of the atmospher es of these planets in transmission 
( iBean et al.ll2011t ICrossfield et al.ll201ll : ICroU et al.ll201ll ). For these observations, precision 
at the level of 1 part in 10^ or better is the objective. At this level, second-order extinc- 
tion effects stemming from variations in molecular absorption may be a dominant source of 
noise, particularly for wide NIR photometric bands. Given a mid-M target and an ensemble 
of G-type reference stars, the theoretical (dashed) curve in the bottom panel of Figure [7] 
indicates the dependence of this effect on PWV. We simulated the impact of PWV varia- 
tions on a long-term M dwarf survey by calculating the differential flux for a mid-M star 
and a mid-G comparison star given the overall distribution of PWV at APO corresponding 
to the times of 6,200 SDSS M star observations. The results of this calculation, shown in 
Figur e [9] indicate that for the SDSS z band and a wide i+z band (INutzman fc Charbonneau 
20081 ). the expected size of this effect is 0.003 to 0.004 mag, comparable to the signal from a 
super-Earth transiting a mid-M star. Large changes in PWV can occur on rapid timescales, 
comparable to the duration of a transit event, further complicating matters. A multiple 
linear regression of the differential flux of an M star target relative to the ensemble reference 
against PWV and AM may reveal significant correlation. The coefficients of this correlation 
are uniquely determined by the SEDs of the target and the reference stars and could be used 
as an empirical relation to improve photornetric measurements by correcting for the impact 
of variations in H2O absorption (see lBurkell201ll ). 



Future wide- area photometric surveys like the Large Synoptic Survey Telescope (LSST; 
Ivezic et al.ll2008l ) will use sophisti cated strateg i es for calibrating the absolute transmission 
of Earth's atmosphere in real time. iBurke et al.l (120 lOf ) demonstrate that spectrophotometric 
observations of a grid of standard stars can be used to calibrate band-integrated atmospheric 
transmission to much better than 1% by modeling both scattering and absorption. While 
absorption by H2O is just one component of the complex and variable atmospheric trans- 
mission function, our results indicate that GPS-based PWV measurements may augment 
spectrophotometric o bservations of standard stars, narrow-band photometric observations 
at key wavelength s (lAlekseeva et al.l l2010al ). or an overall internal calibration procedure 
(j Jones et al.ll2010l ) and help improve photometric precision for objects with strong spectral 
features in the NIR. Another approach to mitigating the effect of variable H2O absorption is 
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to design filters that avoid prominent absorption bands, such as those between 900 and 1000 
nm, though at the cost of the total number of stellar photons detected. Some poss ible designs 
for th e LSST y band filter incorporate a blue edge that avoids H2O absorption (IHigh et al. 
2OIOI). We simulated the impact of variable PWV on differential measurements of a mid-M 
star and a mid-G comparison star using the ys filtei0 shown in Figure [H The results of this 
simulation, shown in Figure [H indicate that careful selection of NIR filters can help suppress 
the impact of variable PWV on M star photometry to below 1 mmag. 

Careful modeling of telluric lines is criti cally important for deriving precise stellar radial 
velocities from high-resolution NIR spectra. IWang &: Gd ( 120111 ) showed that to approach the 
theoretical photon limit in Doppler measurements of low-mass stars, telluric lines need to be 
modeled with residual depths < 1%. While several methods exist for empirical estimation 
of telluric models using observations of A stars (e.g. IVacca et al.ll2003[ ). our results indicate 



that theoretical telluric templates combined with GPS-based PWV data can also be used 
for this purpose. Telluric templates that match observed line depths to better than 1% may 
be an important component of the analysis of data from a new generation of high-resolution 
NIR Doppler instruments (e.g. iMahadevan et al.l 120101 ). 



8. Conclusion 



We calculate theoretical telluric transmission templates in the deep red appropriate for 
Apache Point Observatory (APO) and demonstrate that these models, which contain almost 
exclusively absorption by H2O, are excellent fits to high-resolution observations of A stars. 



These resu lts are quahtat i yely s imilar to those from lThomas-Osip et al.l (120071 ). iQuerel et al. 



( I2OO8I ) and lSeifahrt et al.l ( I2OIOI ) who fit theoretical models to observations of a wide range of 
molecular absorption features in the optical and IR. We compare the scaling parameter for the 
depths of the H2O absorption features around 980 nm, r, to contemporaneous estimates of 
Precipitable Water Vapor (PWV) derived from a Global Positioning System (GPS) receiver 
located at APO and find that they are strongly correlated. A multiple linear regression of 
PWV and airmass (AM) against r indicates that the relative optical depth of H2O can be 
predicted to better than 10%, leading to band-integrated estimates of total water absorption 
that are better than 0.2%. For regions where absorption features have modest optical depth 
and are well-separated, the agreement between the models and telluric absorption features 
in the A star spectra is better than 2% even in the cores of lines. For regions with dense 
forests of hues, leading to saturation, residuals exist at the level of 5 — 10% within the cores 



^ i\protec"t | 



}ittp://ssg. astro. Washington. edu/elsst/magsfilters.shtml?filterinfo 
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of some lines. 

The real-time characterization of atmospheric molecular absorption is important for 
the calibration of broadband NIR photometry. Even using sophisticated differential global 
calibration techniques, changes in PWV can bias measurements of the fluxes of objects 
having strong spectral features in the NIR given that the average SED of the reference 
ensemble of stellar objects is likely to be much smoother. We show that globally calibrated 
SDSS measurements of the z-band fluxes and r-z colors of mid-M stars are correlated with 
PWV. This is not the case for G and F stars, since their NIR SEDs are similar to that 
of the "average" star included in the Ubercal photometric calibration. We compare the 
observed relation between PWV and the M dwarf photometric measurements to theoretical 
predictions based on transmission of the SDSS filters, synthetic stellar templates, and our 
telluric models and find that they are in excellent agreement. We measure the size of the 
impact of changing PWV on SDSS z-band measurements of mid-M stars to be AZH20 = 
(0.0012 ± 0.00008) X [PWV — 2.6 mm) mag. We applied this empirical correction to repeat 
observations from Stripe 82, but found no statistically significant reduction in the variance 
of 25,000 relative z-band measurements of mid-M stars. Given the range of PWV during 
the Stripe 82 observations, the expected size of the correction is (Azhzo) = 0.003 mag, 
significantly smaller than the photon errors for the majority of stars in our sample. 

GPS-based estimates of H2O absorption may be an important component of the calibra- 
tion strategy for fu ture surveys like LSST and, when combin ed with atmospheric monitoring 
( Burke et al.lbolO ) or internal calibration ( Jones et al.lboiol) and a strat e gy for precisel y 



mea- 



suring the spectral response of the instrument (jStubbs fc Tonryl l2006l . IStubbs et al.ll2010l ). 
may help to improve differential photometric measurements for objects with complex NIR 
SEDs. In the near term, modeling H2O absorption through PWV measurements may help 
to improve the precision of differential NIR photometry of bright M stars, enabling transit- 
ing planet surveys to improve their sensitivity to small planets. At the same time, precise 
models of telluric absorption may enable Doppler measurements in new wavelength regimes, 
enhancing ongoing efforts to find Earth-like planets orbiting low-mass stars. 
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Fig. 1. — Components of a tmospheric absorption for a site like Apache Point Observatory 
(see also iBurke et al.ll2010[ ). The transmission curves for NIR bands are also shown, each 
normalized to a peak tr ansmission of 50% . The y band filter shown here is based on the 
proposed LSST ys filter JHigh et al.ll2ninl ). 
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Fig. 2. — Two examples of theoretical telluric H2O templates (solid lines) fit to ARCES A 
star observations (small points). These spectra highlight the strong impact changes in the 
total H2O column have on the telluric absorption. In both cases, the residuals of the fits are 
excellent, having scatter < 2%. 
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Fig. 3. — Two examples of theoretical telluric H2O templates (solid lines) fit to ARCES A 
star observations (small points) in a spectral region with strong H2O absorption, including 
saturated lines. Here we see evidence for significant residuals at the cores of some lines, as 
well as slowly varying residuals that we attribute in part to incomplete correction for fringing 
in the CCD detector. 
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Fig. 4. — Results of a multiple linear regression of line depth scale factor, r, against Precip- 
itable Water Vapor and Airmass. 
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Fig. 5. — Histogram of PWV for ARCES A star observations (shaded region) compared to 
the overall distribution of PWV measurements at Apache Point Observatory. Astronomical 
observations are typically made from APO when PWV< 10 mm. 
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Fig. 6. — Top panel shows mean colors of 6,177 mid-M stars binned as a function of PWV. 
Here, A(r — z) is the difference in the measured color of a star relative to the stellar locus 
in (r-i) vs. (r-z) space. The dashed line is the result of a theoretical prediction that takes 
into account telluric H2O absorptio n scaled following Equation 8, the SDSS filter passbands, 
and synthetic stellar spectra from iBrott fc Hauschildtl (120051 ) . The theoretical (r-z) color 
difference is calculated between a mid-M star and a mid-G star reference, which should be 
representative of the average SED of all the bright stars in SDSS. Error bars are the formal 
statistical error on the mean of each bin and the null hypothesis of A (r — 2;) = is ruled out 
at > 99% confidence. Bottom panel shows results of a consistency check were an empirical 
A{r — z) vs. PWV correction is calculated from a random half of the mid-M stars and then 
this correction is applied to the other half of the observations. After correction, the data are 
consistent with A(r - z) = (P(x^) = 0.73) 
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Fig. 7. — Mean differential z-band magnitude, Az, binned as a function of PWV for objects 
in Stripe 82 that were repeatedly observed by SDSS. The dashed lines are based on the same 
theoretical calculations shown in Figure El again assuming a mid-G star as the reference 
SED. Formal statistical errors on the means are smaller than the points. 14,597 observations 
of 983 stars are included in the top panel and 25,005 observations of 1,553 stars in the bottom 
panel. 
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Fig. 8. — Mean values of 14,597 uncalibrated z-band flux measurements for the Stripe 82 G5- 
F5 star sample binned as a function of PWV. As discussed in Section 6, ^ Zraw is an estimate 
of the change in raw z-band flux after correction for extinction as measured simultaneously 
in r band. The dashed line is based on the same theoretical calculations shown in Figure El 
The changes in total absorption in z-band with PWV are significant, up to several percent, 
but these are largely removed via relative calibration techniques, with residuals that depend 
on SED. Statistical errors on the means are smaller than the points. 
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Fig. 9. — Histograms of the expected variation in differential magnitude for a mid-M target 
and mid-G comparison star in different NIR bands given tlie PWV values at the times of all of 
the SDSS observations in our M star sample. The standard deviations of these distributions 
are: aAz = 0.003, aAi+z = 0.004, and a^y = 0.0003 mag. 
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Table 1. Empirical relation between PWV measured at White Sands and PWV measured 

at APO. 
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